Asymptotic regime in N random interacting species 
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cn ■ Abstract 



The asymptotic regime of a complex ecosystem with N random interacting species and in the 
presence of an external multiplicative noise is analyzed. We find the role of the external noise on 



presence or the absence of some population give different asymptotic distributions of these fields. 
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I. INTRODUCTION 



In recent years great attention has been devoted to population dynamics modelled by 
generalized Lotka-Volterra systems JJ. Ecosystems are a classic example of complex sys- 
tems, which became object of study as well by biologists as by physicists P, [|. Tools 
developed in the context of nonequilibrium statistical physics to analyze nonequilibrium 
nonlinear physical systems provide new insights and at the same time new approaches to 
the comprehension of the properties of biological and many body systems. A key aspect 
to understand the complex behavior of ecosystems is the role of the external noise on the 
dynamics of such systems. Noise-induced effects in population dynamics, such as pattern 
formation , stochastic resonance, noise delayed extinction, quasi periodic oscillations 
etc... have been investigated with increasing interest 

mm 

. The dynamical behavior of 

ecological systems of interacting species evolves towards the equilibrium states through the 
slow process of nonlinear relaxation, which is strongly dependent on the random interaction 
between the species, the initial conditions and the random interaction with environment. 
Moreover biological evolution presents the same fundamental ingredient that characterizes 
physical systems far from equilibrium in their route to equilibrium, namely the disorder- 
order transition. Different models of evolution are reported in literature, which are useful 



to describe a lot of dynamical population problems 
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llj . Among them it is worthwhile 
to cite two bit-string models of population dynamics, namely the Eigen model [12] and the 
Penna model jl^ ]. 

The mathematical model here used to analyze the dynamics of N biological species with 
spatially homogeneous densities is the generalized Lotka-Volterra system with a Malthus- 
Verhulst modelling of the self regulation mechanism and with the addition of an external 
multiplicative noise soorce flfi. We obtain the a 8y ,nptot 1C behaviors of the probability 
distribution of the populations for different values of external noise intensity. To analyze 
the role of the external noise on the transient dynamics of the species we focus on the long 
time distribution of the local field, that is the interaction term in the dynamical equation 
of the i th species that takes the influence of all other species into account. We find that 
the presence or the absence of some population give different asymptotic distributions of 
the local field and of the cavity field (field acting on the i th species when this is absent) 
in the absence of external noise. When the noise is switched on the asymptotic local and 
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cavity fields tend to overlap and approximately superimpose each other for very high noise 
intensity. Finally the long time evolution of the average number of the extinct species is 
reported for different values of the multiplicative noise intensity. 



II. THE MODEL 



The dynamical evolution of our ecosystem composed by TV interacting species in a noisy 
environment (climate, disease, etc..) is described by the following generalized Lotka-Volterra 
equations with a multiplicative noise, in the framework of Ito stochastic calculus 



drii (t) 



IWj 



7lj(t), 



i = l,...,N 



1) 



gi(ni(t)) + ^ I dt + Vtd 

where nj(t) > is the population density of the i th species at time t and the function gi(rii(t)) 



g i (n i (t))=(a + ^-n i (t), (2) 

describes the development of the i th species without interacting with other species. In 
Eq. (JJJ, a is the growth parameter, the interaction matrix models the interaction between 
different species (i ^ j), and Wi is the Wiener process whose increment dwi satisfy the usual 
statistical properties (dwi(t)) = 0, and (dwi(t)dwj(t')) = 5ij5(t — t')dt. We consider an 
asymmetric interaction matrix J^, whose elements are randomly distributed according to 
a Gaussian distribution with (Jij) = 0, (JijJji) = 0, and <r| = J 2 /N. Therefore our 
ecosystem contains 50% of prey-predator interactions ( < and Jjt > 0), 25% competitive 
interactions (Jy > and Jji > 0), and 25% symbiotic interactions (Jy < and Jji < 0). 
We consider all species equivalent so that the characteristic parameters of the ecosystem are 
independent of the species. The formal solution of Eq. (JIJ) is 



Ui{t) 



ni(0)zi{t) 



1 + ^(0) f*dt' Zi (t') 



where 



i(t) = exp at + y/ewi(t) + / dt'h iMc (t') 

Jo 



(3) 



(4) 



The term hi t i oc {t) = Ylj^i Jij n j(f) is the local field acting on the i th population and represents 
the influence of other species on the differential growth rate. We note that the dynamical 
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behavior of the i th population depends on the time integral of the process Zi(t) and the time 
integral of the local field. 

In the absence of external noise (e = 0), for a large number of interacting species we can 
assume that the local field hi(t) is Gaussian with zero mean and variance 

al iioc = T ljjk (J ij J ik )(n j n k ) = J 2 {nj) , with (JijJik) = 5 jk — . (5) 

As a consequence, in the absence of external noise, from the fixed-point equation rii(a — 
ni + hi) = 0, the stationary probability distribution of the populations is the sum of a 
truncated Gaussian distribution at n« = (jii > always) and a delta function for extinct 
species. The initial values of the populations rij(0) have also Gaussian distribution with 
mean value (rij(0)) = 1, and variance &n(o) = 0-03. The interaction strength J determines 
two different dynamical behaviors of the ecosystem. Above a critical value J c = 1.1, the 
system is unstable and at least one of the populations diverges. Below J c the system is stable 
and asymptotically reaches an equilibrium state. The equilibrium values of the populations 
depend both on their initial values and on the interaction matrix. If we consider a quenched 
random interaction matrix, the ecosystem has a great number of equilibrium configurations, 
each one with its attraction basin. For an interaction strength J = 1 and an intrinsic growth 
parameter a = 1 we obtain: (rii) = 1.4387, (nf) = 4.514, and o~\, = 2.44. These values agree 
with that obtained from numerical simulation of Eq. (Q). 

In the presence of external noise (e 7^ 0) we calculate long time probability distribution for 
different values of the noise intensity. These are shown in the following Fig. ^ For increasing 
external noise intensity we obtain a larger probability distribution with a lower maximum 
(see the different scales in Figs. [0 for different noise intensity values). The distribution 
becomes asymmetric for e = 0.1 and tends to become a truncated delta function around the 
zero value (-P(^j) = 5{n{) for rij > 0, and P(rii) = for rij < 0), for further increasing noise 
intensity. The role of the multiplicative noise is to act as an absorbing barrier at = ^| . 
To confirm this picture we calculate the time evolution of the average number of extinct 
species for different noise intensities. This time behavior is shown in Fig. |21 We see that 
this number increases with noise intensity, and after the value e = 0.1 reaches quickly the 
normalized maximum value at e = 10. 

To analyze in more detail the influence of each species on the long time dynamics of 
the ecosystem we calculate in this regime the local field h it i oc (t) and the cavity field h itCav (t), 
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FIG. 1: Probability distribution for the species densities. The values of the external noise intensity 
are e = 0,0.001,0.01,0.1, 1, 10. 
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Number of Extinct Species: N ext (t) 
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FIG. 2: Time evolution of the average number of extinct species for different noise intensities. 



which is the field acting on the i th population when this population is absent. The probability 
distributions for both local and cavity fields have been obtained by simulations for different 
species in the presence and in absence of external noise. The results are shown in the next 
section (see Fig. |3J). We found that the probability distributions of the cavity fields differ 



5 



substantially from those of local fields for the same species, while in the presence of noise 
the two fields overlap. To quantify this overlap between the probabilities distributions of the 
two fields we define an overlap coefficient A(t), which is the distance between the average 
values of the two distributions, normalized to their widths 

\/,\ hii oc hi cav dhit) , , 

2 + 2 a d {t) 

iloc 1 i.cav 



where 

d h (t) = hi,loc ~ hcav , °l{t) = vlloM + vlcavit) ■ ( 7 ) 

With this definition the distributions start to overlap significantly for |A| < 1, and become 
strongly overlapping for |A| <C 1. 



III. RESULTS AND COMMENTS 

In the calculation the following parameters have been used: a = 1.2, J = 1, ctj = 
0.005, N = 200; the number of averaging experiment used is N exp = 1000. Concerning 
the initial condition the parameters are: (n^) = 1, <r^ o = 0.03. The dynamics of various 
species are different even if they are equivalent according to the parameters in the dynamical 
Eq. (^) . However we note that to change the species index by fixing the random matrix or to 
change the random matrix by fixing the species index is equivalent as regards the asymptotic 
dynamical regime. Fig. El shows for various noise intensities the local and cavity fields 
probability distributions at time t = 50 (a. u.). For noiseless dynamics the distributions of 
the fields for the species % = 3 appear very narrow around their mean values and very spaced 
each other. By increasing the noise intensity, we observe an equal enlargement of the two 
distributions, maintaining however the same mean values. At e = 0.1 the two distributions 
start to overlap until, for stronger noise intensity (e = 1), they superimpose each other. The 
overlap coefficient \(t) is equal to zero. The noise makes equivalent all the species in the 
asymptotic regime and the absence of some species doesn't contribute to any changes in the 
dynamics of all other species. The last plot in Fig. El (e = 10) gives a delta distribution 
around zero. This means that, at the time considered in our simulation (t = 50) and for 
this noise intensity, all the species are extinct (see also Fig- EJ • 

The detailed time evolution of the overlap of the two distributions can be seen from 
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FIG. 3: Probability distribution of both the local (black circles)vand the cavity (white circles) 
fields for various values of noise intensity e = 0, 0.001, 0.01, 0.1, 1. The graph are taken at the time 
t = 50 unit steps. 



Fig. EJ where it is plotted the coefficient denned in Eq. (JHJ). For e = the X(t) coefficient 
increases with time. This is due to the different time behavior of the distance between the 
mean values of the field distributions and of their standard deviation af loc and crf cav . The 
distance dh{t) = hi t i oc — hi >cav is almost constant in time, except a rapid initial transient (see 
Fig- El e = 0), but the corresponding evolution of the distribution widths decreases rapidly in 
time. This effect is due to the quenched random matrix. This behavior remains unchanged 
until the noise intensity reaches the value of e = 0.01. ^From this value of external noise 
intensity some differences start to be visible (see Figs. HJ E] and EI), and at e = 1, after some 
fluctuations, both the distance dh(t) and the overlap coefficient X(t) reach a value close to 
zero. The two field distributions are totally overlapped. 

It is worthwhile to note that the behaviors shown in Figs. El lU El and El have been found 
for some species, and changing the species different evolutions of the distribution dynamics 
appear and in particular of the A(t) coefficient. This is due to complexity of our ecosystem 
and to the extinction process during the transient dynamics. Moreover this strange behav- 
ior, found for some populations and in the asymptotic regime, is reminiscent of the phase 
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FIG. 4: Time evolution of the overlap coefficient X(t) between the distributions of local and cavity 
fields for different noise intensities. The increasing behavior of the overlap coefficient as a function 
of time, in the noiseless case, is due to the narrowing of the two fields distributions (o7 oc , a cav — ► 0) 
towards a 5-function. For a noisy environment the distributions tend to enlarge, decreasing the 
value of the overlap coefficient X(t). 

transition phenomenon [l^j], and it is related to the following peculiarities of our dynamical 
system: (i) all the populations are positive; (ii) different initial conditions drive the ecosys- 
tem into different attraction basins; and (iii) the complex structure of the attraction basins. 
While in the presence of noise all the populations seem to be equivalent in the long time 
regime, some populations, in the absence of external noise, have an asymptotical dynami- 
cal behavior such that they significantly influence the dynamics of other species. A more 
detailed analysis on these points will be done in a forthcoming paper. 



IV. CONCLUSIONS 



We analyzed the asymptotic regime of an ecosystem composed by N interacting species 
in the presence of multiplicative noise. We find the role of the noise on the asymptotic 
probability distribution of populations and on the extinction process. Concerning the local 
and the cavity fields, a phase transition like phenomenon is observed. Their probability 
distributions tend to overlap each other in the presence of external noise, reaching strong 
overlap for high noise intensity (|A(t)| ~ 0), while they are separated (|A(t)| > 1) in the 
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FIG. 5: Time evolution of the distance dh(t) between the mean values of the fields distributions 
for different noise intensities, namely e = 0,0.001,0.01,0.1, 1. 
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FIG. 6: Time evolution of the standard deviation crf loc (grey circles)and <jf cav (black circles) of 
the field distributions for the same noise intensities of Fig. [21 



absence of noise. This phenomenon can be ascribed to the complexity of our ecosystem. 



9 



Acknowledgements 

This work was supported by MIUR and INFM-CNR. 



[1] G. Baxter, A. J. McKane, and M. B. Tarlie, Phys Rev. E 71, 011106 (2005); A. J. McKane and 
T. J. Newman, Phys. Rev. Lett. 94, 218102 (2005); A. Shabunin, A. Efimov, GA. Tsekouras, 
et al., Physica A 347, 117-136 (2005); K. Tokita, Phys. Rev. Lett. 93, 178102 (2004); G. J. 
Ackland and I. D. Gallagher, Phys. Rev. Lett. 93, 158701 (2004); R. Mankin, A. Sauga, A. 
Ainsaar, et al., Phys. Rev. E 69, 061106 (2004); Y. De Decker, GA. Tsekouras, A. Provata, 
et al, Phys. Rev. E 69, 036203 (2004); GA. Tsekouras, A. Provata, C. Tsallis, Phys. Rev. 
E 69, 016120 (2004); A. Provata, GA. Tsekouras, F. Diakonos, et al., Fluct. Noise Lett. 3, 
L241-L250 (2003). 

[2] J. D. Murray, Mathematical Biology, 2nd ed. (Springer, New York, 1993) 54 (1996). 

[3] C. Escudero, J. Buceta, F. J. de la Rubia, and Katja Lindenberg, Phys. Rev. E 69, 021908 

(2004); T. J. Kawecki and R. D. Holt, Am. Nat. 160, 333 (2002); Michel Droz and Andrzej 

Pekalski, Phys. Rev. E 69, 051912 (2004); I. Giardina, J. P. Bouchaud, M. Mezard, J. Phys. 

A: Math. Gen. 34, L245 (2001); H Rieger, J. Phys. A: Math. Gen. 22, 3447 (1989). 
[4] J. E. S. Socolor, S. Richards, and W. G. Wilson, Phys. Rev. E 63, 041908 (2001). 
[5] A. Fiasconaro, D. Valenti and B. Spagnolo, Acta Phys. Pol. B 35, 1491 (2004); D. Valenti, 

A. Fiasconaro and B. Spagnolo, Acta Phys. Pol. B 35, 1481 (2004); A. La Barbera and B. 

Spagnolo, Physica A 314, 120 (2001). 
[6] B. Spagnolo A. Fiasconaro and D. Valenti, Fluct. Noise Lett. 3, L177 (2003); B. Spagnolo 

and A. La Barbera, Physica A 315, 201 (2002); A. F. Rozenfeld Rozenfeld, C.J. Tessone, E. 

Albano, H.S. Wio, Phys. Lett. A 280, 45 (2001); J. M. G. Vilar and R. V. Sole, Phys. Rev. 

Lett. 80, 4099 (1998). 

[7] See the special section on "Complex Systems", Science 284, 79-107 (1999); the special section 

on "Ecology through Time", Science 293, 623-657 (2001). 
[8] D. Valenti, A. Fiasconaro and B. Spagnolo, Physica A 331, 477 (2004). 

[9] B. Spagnolo D. Valenti, A. Fiasconaro, Math. Biosciences and Eng. 1, 185 (2004); D. F. 
Russel, L. AQ. Wilkens and F. Moss, Nature 402, 291 (2000). 



10 



[10] S. Moss de Oliveira, Domingos Alves, and J. S. Sd Martins, Physica A 285, 77 (2000). 
[11] LucaPeliti, "Introduction to the statistical theory of Darwinian evolution", cond-mat/9712027, 
1997. 

[12] M. Eigen, Naturwissenschaften, 58, 465 (1971); M. Eigen, J. Mc Caskill, and P. Schuster, 

Adv. Chem. Phys. 75, 149 (1989). 
[13] T. J. P. Penna, J. Stat. Phys. 78, 1629 (1995). 

[14] S. Chichi, F. de Pasquale and B. Spagnolo, Phys. Rev. E 54, 706 (1996); ibid. 47, 3915 (1993); 

P. Barrera, S. Ciuchi and B. Spagnolo, J. Phys. A: Math. Gen. 26, L559-L565 (1993). 
[15] B. Spagnolo, M. A. Cirone, A. La Barbera and F. de Pasquale, J. Phys.: Condens. Matter 

14, 2247 (2002); M. A. Cirone, F. de Pasquale and B. Spagnolo, Fractals 11, 217 (2003); B. 

Spagnolo, D. Valenti and A. Fiasconaro, Prog. Theor. Phys. Suppl. 157, 312-316 (2005). 
[16] S. Ciuchi, F. de Pasquale, P. Monachesi and B. Spagnolo, Phys. Rev. B 38, 2596 (1988); S. 

Ciuchi, F. de Pasquale, P. Monachesi and B. Spagnolo, Physica Scripta vol. T25, 156 (1989). 



11 



